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ABSTRACT 

Context. A significant fraction of cosmological gamma-ray bursts (GRBs) are characterised by a fast rise and exponential decay 
(FRED) temporal structure. This is not a distinctive feature of this class, since it is observed in many Galactic transients a nd is likely 
descrip tive of a sudden release of energy followed by a diffusion process. Possible evidence has recently been reported bv lTello et alj 
ll2012h for a Galactic contamination in the sample of FRED GRBs discovered with Swift. 

Aims. We searched for possible Galactic intruders disguised as FRED GRBs in the Swift catalogue up to September 2014. 

Methods. We selected 181 FRED GRBs (2/3 with unknown redshift) and considered different subsamples. We tested the degree of 
isotropy through the dipole and the quadrupole moment distributions, both with reference to the Galaxy and in a coordinate-system— 
independent way, as well as with the two-point angular autocovariance function. In addition, we searched for possible indicators of a 
Galactic origin among the spectral and temporal properties of individual GRBs. 

Results. We found marginal (~3cr) evidence for an excess of FREDs with unknown redshift towards the Galactic plane compared 
with what is expected for an isotropic distribution corrected for the non-uniform sky exposure. However, when we account for the 
observational bias against optical follow-up observations of low-Galactic latitude GRBs, the evidence for anisotropy decreases to 
~2(T. In addition, we found no statistical evidence for different spectral or temporal properties from the bulk of cosmological GRBs. 
Conclusions. We found marginal evidence for the presence of a disguised Galactic population among Swift GRBs with unknown 
redshift. The estimated fraction is / = (19 ± 11)%, with an upper limit of 34% (90% confidence). 

Key words, gamma-ray burst: general - methods: statistical 


1. Introduction 


Gamma-ray bursts (GRBs) are known to be cosmological tran¬ 
sient sources that signal the very final stage of some kind of 
massive stars for the long duration ones, and likely the act of 
mer ging for a compact binary system for short duration ones 
fsee lMeszaros & Gehrelsll2012l for a review). Owing to their ex¬ 
treme luminosity they are routinely detected from cosmol ogical 
distances up to redshifts z ~ 9.4 (ICucchiara et al.ll201 ll) . The 
consequent highly isotropic sky distribution was interpreted as 
evidence f or a cosmological origin prior to the first redshift mea¬ 
surements (iHartmann & Epsteinlll9^lHartmann & Blumenth^ 
1989HBriggs et al.ll996l) . Once the cosmological nature of long- 

duration GRBs was finally established, a number of studies 
reported eviden ce for anisotropy i n the distribution of short- 
duration ones (iBalazs et al.1119981: Magliocchetti et al.l 120031: 
IVavrek et al.ll2008 ; but see also iBernui et al.l ^08h . or evidence 
for a corr elation between sh ort GRBs and galaxies in the local 
Universe (iTanvir et al.ll2005l) . 

Recently, Tello et al. (2012; hereafter T12) exploited 
the real-time arcmin-sized lo calisation capabilities o f the Swift 
Burst Alert Telescope (BAT; iBarthelmv et al.l 120051) to test the 
isotropy of a special class of long GRBs, the so-called fast-rise 
exponential-decay (ERE D) ones. Not only i s their time profile 
quite common in GRBs (iNorris et al.lll996l) . but it is also de¬ 
scriptive of various high -energy outbursts from Galactic sources, 
such as X-ray binaries dRemillard & McClintockll2006h . X-ray 
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bursters (iLewin et akl 119931) . magnetars (iMereghe ttil 2008), as 
well as unidentified Galactic trmsients d ^asliwal et al.l 20081: 
ICastro-Tirado et aLll2008l: IStefanescu et al, l2008l) . Such univer- 
sality is observed across different wavelengths all the way to the 
radio bands, where analogous flares are observed over timescales 
from seconds or les s, up to years, due to a broad range of as- 
trophysical sources (iPietka et al.11^1^ . The ubiquity of ERED- 
like transient profiles is explained by its being descriptive of a 
sudden energy release, followed by some cooling and diffusion 
process. T12 found some evidence (< 3cr confidence) for a devi¬ 
ation from isotropy of the sky distribution of Swift EREDs with 
unknown z. They attribute it to the presence of Galactic uniden¬ 
tified sources which would make up to ~ 27% of the observed 
unknown-z ERED population. 


Motivated by their results, we carried out a similar investiga¬ 
tion on an updated sample of Swift ERED catalogue, adopting in¬ 
dependent selection criteria. We analysed the degree of isotropy 
both with reference to the Galactic plane and in a reference- 
system-independent way. In addition, to further characterise the 
nature of the Galactic candidates, we examined on a statistical 
basis the spectral and temporal properties of the y-ray prompt 
and X-ray afterglow e mission as observed with Swift BAT and 
X-ray Telescope (XRT: [Burrows et al.l2005l) and compared with 
a sample of cosmological GRBs. The paper is organised as fol¬ 
lows: the sample selection and data analysis are described in Sec- 
tions|2]and[2 respectively. Section|4]presents the results, which 
are discussed in Section [5] 


Article number, page 1 of|6] 




















































:’red_v3 



00 

o 

o 

CD Fig- 1- Exposure map of BAT in Galactic 
o coordinates from January 2005 to Septem¬ 
ber 2014. Green diamonds (blue circles) are 
Q all FREDs with unknown (known) z. They 
jointly constitute SI. All green diamonds 
are sample S2. Crossed green diamonds are 
those unknown-z FREDs with a unique de- 
cD tected peak (sample S3). The three sets in- 
o elude 181, 119, and 71 GRBs, respectively. 


Table 1. GRB samples’ sizes. 


Sample 

Description 

Size 
(this work) 

Size 

(T12) 

SI 

all FREDs 

181 

111 

S2 

all FREDs, no z 

119 

77 

S3 

all FREDs, no z, 1 peak 

71 

49 

STot 

all GRBs 

904 

- 

SnFnZ 

all non-FREDs, no z 

506 

- 


2. Sample selection 

Starting from a sample of 904 GRBs detected by 5w!/f-BAT in 
the time interval January 2005 to September 2014, we selected 
those which had explicitly been tagged by the BAT team as 
FREDs in the GCN circulars. We ended up with 181 FREDs out 
of 904 GRBs, whereas T12 found 111 FREDs out of 596 GRBs. 
We further characterised the number of pea ks of each GRB by 
means of the MEPSA code (lGuidorzill2015h . capable of identi¬ 
fying relatively faint peaks at very different timescales. We set a 
threshold of 4.5 on the minimum signal-to-noise as yielded by 
MEPSA, which ensures a negligible false positive rate for our 
light curve sample. We then considered three different samples: 
i) the full one; ii) the subsample of GRBs with unknown red- 
shift z; iii) the subsample of GRBs with unknown z and with a 
single peak detected with MEPSA. Hereafter these samples are 
referred to as SI, S2, and S3, including 181, 119, and 71 GRBs, 
respectively. Their sky distributions in Galactic coordinates are 
shown in Fig. [T] To assess the role of observational biases, for 
comparison we also consider two additional samples detected in 
the same time interval: (i) all 904 BAT-discovered GRBs, and 
(ii) the complementary sample of 506 non-FREDs GRBs with 
unknown z. Hereafter, we refer to them as STot and SnFnZ, re¬ 
spectively. Table[T]summarises the samples’ sizes and the corre¬ 
sponding values of T12. 

3. Data analysis 

For each real sample we calculated the three Galactic dipole 
moments (sin b, cos b sin /, cos b cos 1), and the five quadrupole 
moments (sin^Z? - 1/3, cos^Z7sin2/, cos^bcosll, sin2fesin(, 
sin 2b cos /) with an emphasis on the ones related to the Galactic 
plane, i.e. sinZ? and (sin^ b - \ /3). We also calculated the some¬ 
what redundant equatorial dipole sin 6 and quadrupole (sin^ 6 - 
1 /3) moments to better evaluate possible anisotropies connected 
with observational biases. In addition , we calculated the two- 
point angular correlation function dHartmann & Blumenthall 


Il989l) . We also carried out a reference-system-independent 
search for anisotropies possibly unrelated with our Galaxy by 
eval uating the Ra yleigh-Watson 'W and the Bingham S statis¬ 
tics (lBriggslll993[) . Everything was then compared with the ex¬ 
pectations for isotropy following two independent approaches: 
through Monte Carlo simulations and sky pixelisation. 

The exposure map of 5w!/f-BAT (Fig. [T]| was obtained from 
the limiting flux map for the considered time inte rval, process¬ 
ing th e BAT data with the Bat_lmager softwareJSegreto_et^ 
I 2 OIOI). and following the procedure described in ICusumano et al. 

3.1. Monte Carlo simulations (method 1) 

For each sample we generated 10^ synthetic samples, each hav¬ 
ing as many positions as the corresponding real one. Each fake 
position was generated according to the following procedure: i) 
a random position is drawn from an isotropic distribution; ii) 
the position is accepted/rejected depending on the outcome of 
a binomial variate, whose acceptance probability was set to the 
fractional exposure of that position. The procedure ends as soon 
as the desired number of positions is achieved. 

3.2. Sky pixelisation (method 2) 

The entire sky was split into a homogen eous grid of 529 2 pix¬ 
els following the technique described by iTegmarkI (Il996h . with 
angular resolution of 3?2. Each pixel was then assigned the cor¬ 
responding fractional exposure value. The choice of the angular 
resolution was driven by the need to ensure an adequate cov¬ 
erage of the BAT exposure map. The expected moments were 
calculated over the entire sky grid by weighting the contribution 
of each sky pixel by the associated fractional exposure. The un¬ 
certainties affecting the observed moments were calculated from 
the expected variances in the Gaussian limit ensured by the cen¬ 
tral limit theorem due to A » 1 events: cr^((sinf7)) = 1/(3A), 
o-2«sin^ b)) = 4/(45A) for the ( sin^ -1/3) t erm, cr^ = 4/(15A) 
for the other quadrupole terms (lBriggslll993l) . 

4. Results 

4. 1. Dipole and quadrupole moments 

Table |2] reports the observed moments along with the expected 
values and intervals. While none of the dipole moments deviate 
by >2cr from isotropy, the Galactic quadrupole moment of S2 is 
marginally (-2.8cr) lower than the isotropic-expected one. Anal¬ 
ogous result is obtained for SnFnZ with a less significant devi¬ 
ation though (-2.1cr). This suggests that the unknown-z FREDs 
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Table 2. Average dipolar and quadrupolar moments for the ditferent samples. The real value moments are reported and associated to the values 
obtained from the sky pixelisation method, along with their discrepancies in cr units. The last two columns report the MC Icr and 2cr (Gaussian- 
equivalent) confidence intervals. 


Sample 

Moment 

Observed 

Expected^"-* 
(sky pix.; +lcr) 

Discrepancy 

(o-) 

Expected Interval (MC)^*^ 

(Icr) (2cr) 

SI 

(cos b cos 1) 

-0.010 

-0.010 + 0.043 

0.0 

-0.050, +0.031 

-0.091,+0.072 

SI 

{sinb) 

0.028 

0.022 + 0.043 

+0.1 

-0.018,+0.064 

-0.061,+0.105 

SI 

(sin 6) 

0.065 

0.035 + 0.043 

+0.7 

-0.012, +0.074 

-0.052,+0.120 

SI 

(sin^ b - 1/3) 

-0.009 

0.014 + 0.022 

-1.0 

-0.009, +0.035 

-0.030, +0.057 

SI 

(sin^^- 1/3) 

0.074 

0.039 + 0.022 

+ 1.5 

+0.016, +0.062 

-0.005, +0.085 

S2 

(cos b cos 1) 

-0.005 

-0.010 + 0.053 

+0.1 

-0.063, +0.036 

-0.108,+0.081 

S2 

{sinb) 

0.097 

0.022 + 0.053 

+ 1.4 

-0.030, +0.075 

-0.089,+0.125 

S2 

(sin 6) 

0.070 

0.035 + 0.053 

+0.7 

-0.018,+0.086 

-0.072,+0.141 

S2 

(sin^(7- 1/3) 

-0.063 

0.014 + 0.027 

-2.8 

-0.012, +0.043 

-0.036, +0.072 

S2 

(sin^^- 1/3) 

0.103 

0.039 + 0.027 

+2.3 

+0.011,+0.068 

-0.016,+0.098 


(cos b cos /) 

-0.005 

0.017 + 0.053 

-0.4 

-0.036, +0.066 

-0.080, +0.123 

S2^"'' 

(sinZt) 

0.097 

-0.008 + 0.053 

+2.0 

-0.057, +0.043 

-0.108,+0.096 

S2^"'’ 

(sin 6) 

0.070 

-0.022 + 0.053 

+ 1.7 

-0.079, +0.034 

-0.130, +0.085 

S2^'^'' 

{sm^b-l/3) 

-0.063 

-0.013 + 0.027 

-1.8 

-0.042, +0.015 

-0.069, +0.045 

S2^'^'' 

(sin^^- 1/3) 

0.103 

0.056 + 0.027 

+ 1.7 

+0.031,+0.086 

-0.001,+0.114 

S3 

(cos b cos /) 

-0.005 

-0.010 + 0.068 

+0.1 

-0.080, +0.059 

-0.148,+0.117 

S3 

(sinZ?) 

0.111 

0.022 + 0.068 

+ 1.3 

-0.046, +0.088 

-0.109,+0.165 

S3 

(sin 6) 

0.097 

0.035 + 0.068 

+0.9 

-0.038,+0.106 

-0.109,+0.171 

S3 

{sm^b-l/3) 

-0.032 

0.014 + 0.035 

-1.3 

-0.025, +0.048 

-0.059, +0.083 

S3 

(sin^^- 1/3) 

0.102 

0.039 + 0.035 

+ 1.8 

+0.004, +0.078 

-0.030, +0.121 

SnFnZ 

(cos b cos /) 

0.036 

-0.010 + 0.026 

+ 1.8 

- 

- 

SnFnZ 

(sinZt) 

0.011 

0.022 + 0.026 

0.0 

- 

- 

SnFnZ 

(sin 6) 

0.040 

0.035 + 0.026 

+2.0 

- 

- 

SnFnZ 

{sin^b - 1/3) 

-0.014 

0.014 + 0.013 

-2.1 

- 

- 

SnFnZ 

(sin^^- 1/3) 

0.046 

0.039 + 0.013 

+0.5 

- 

- 

STot 

(cos b cos 1) 

0.014 

-0.010 + 0.019 

+ 1.2 

- 

- 

STot 

(sinb) 

0.011 

0.022 + 0.019 

-0.6 

- 

- 

STot 

(sin 6) 

0.040 

0.035 + 0.019 

+0.2 

- 

- 

STot 

(sin^b - 1/3) 

0.017 

0.014 + 0.010 

+0.4 

- 

- 

STot 

(sin^cJ- 1/3) 

0.037 

0.039 + 0.010 

-0.2 

- 

- 


Notes. Estimated through the sky pixelisation (method 2) in the Gaussian limit ensured by the central limit theorem. Estimated through 
MC simulations (method 1). S2 is compared with expectations that modelled observational biases connected with the redshift measurement 
(Section lA?! . 
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sin^(b) -1/3 


Fig. 2. Galactic quadrupole moment distributions obtained from MC simulations (method 1; histograms) and from the sky pixelisation in the 
Gaussian limit (method 2; solid curves) compared with the values observed in the real samples (vertical lines) for SI, S2, and S3 (left, mid, and 
right panels, respectively). 
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tend to cluster at low Galactic latitudes, or equivalently, that 
they are rarely found at high latitudes. Among the other Galactic 
quadrupole moments, only sin 2b cos / exhibits two > 2cr devia¬ 
tions from isotropy, -2.5cr (-2.4cr) for SI (S2). Figure |2] shows 
the Galactic quadrupole moment (sin^ b - I /3) for each of the 
three samples obtained with both methods assuming isotropy, 
corrected for the non-uniform exposure map, to be compared 
with the value observed in the corresponding real sample. 


4.2. ^ and S statistics 

Following iBriggsl (Il993h . for each position (/,, bi) we consid¬ 
ered the associated unit vector r,- = {Xi,yi,Zi) and calculated the 
Rayleigh-Watson statistic defined as 




( 1 ) 


This statistic characterises the dipole moment and in the absence 
of any preferred direction is asymptotically distributed as for 
large N. The Bingham statistic S is defined as 



where A/; are the eigenvalues of the so-called “orientation” ma¬ 
trix, that is equivalent to the quadrupole one and is defined as 


l=l 


' XiXi 

xtyi 

XiZi ' 

yiXi 

yiyt 

yiZi 

. ZiXi 

Ziyi 

ZiZi > 


(3) 


Table 3. Observed values of Rayleigh-Watson (TV) and Bingham (S) 
statistics and confidence intervals expected for an isotropic distribution 
corrected for the non-uniform exposure. 


Sample Statistics Observed Expected Interval (MC)'* 

(Icr) (2o-) 


SI 

TV 

2.6 

0.9 

6.0 

0.2 

10.5 

S2 

TV 

3.7 

0.9 

5.6 

0.2 

10.2 

S3 

TV 

3.3 

0.9 

5.5 

0.2 

10.1 

SI 

S 

14.2 

4.9 

16.0 

1.7 

24.8 

S2 

S 

19.4 

3.7 

13.8 

1.6 

21.5 

S3 

tB 

9.8 

2.9 

11.1 

1.1 

18.6 


Notes. Estimated through MC simulations (method 1). 



Angular distance 6 [deg] 


S measures the deviation of the eigenvalues from the value 1 /3 
expected for isotropy and is asymptotically distributed as Xy Be¬ 
cause of the non-uniform sky exposure we applied just method 1 
to derive the expected intervals for both statistics. TableOreports 
the observed values and the corresponding intervals for isotropy. 
Like for the Galactic moments we did not find evidence for any 
preferred direction with a non-zero dipole moment, all within 1 cr 
confidence. As for the quadrupole term, although all of them lie 
within 2cr, S2 shows the largest deviation. In particular, it mildly 
suggests a clustering of unknown-z FREDs along the Galactic 
direction / = 130°, h = -hl3° (a = 39°, 5 = h- 75°, J2000) and 
its antipodal one, i.e. close to celestial poles. Erom Table |2] the 
equatorial quadrupole moment (sin^ 5-1 /3) exceeds by 2.3cr the 
expected isotropic value only for S2, thus lending support to the 
enhanced presence of unknown-z EREDs towards the celestial 
poles. 

The unknown-z non-EREDs GRBs (SnEnZ) show a similar 
behaviour: their Galactic quadrupole moment deviates by -2.1cr 
from isotropy. In contrast, in spite of the best statistical sensitiv¬ 
ity the whole BAT-detected sample (STot) shows no evidence at 
all for non-zero quadrupole moments. 


Fig. 3. Histogram of the angular autocovariance function distribution 
observed for S2 compared with what is expected taking into account 
the non-uniform sky exposure (thick red line). The bin width is uniform 
in cos 8. 

most remarkable deviations are observed for S2, as shown in 

Fig.E] 

The thick line shows the isotropic-expected value as a func¬ 
tion of the angular distance. The displayed uncertainties on the 
observed numbers are their square root values, under the as¬ 
sumption of Poisson statistics. Strictly speaking, this is not cor¬ 
rect: the ensemble of 6 values are not fully independent of each 
other, since the entire set of possible pairs is not. However, given 
the large number of pairs the degree of correlation is relatively 
small, so as a matter of fact Poissonian uncertainties are realistic. 

Unknown-z EREDs seem to cluster on angular scales < 30° 
(between 2cr and 3cr significance) at the expense of those in the 
range 60°-90°. A x^ test yields a p-value of 15% {x^ fdof - 
25.4/19). However, this must be taken as a loose indication, for 
the lack of statistical independence of the different bins. 


4.3. Two-point anguiar autocovariance 


4.4. Observationai biases connected with redshift 
measurement 


Einally, we calculated the distribution of the two-point angular 
autocovariance function by taking all the pairs within a given 
sample. In particular, we considered cos 9 rather than 6 (angu¬ 
lar distance between a pair), given that for isotropy and uniform 
exposure the expected distribution is uniform in cos 0. Likewise, 
the expected distribution for isotropy was derived only through 
method 1 to account for the non-uniform exposure map. The 


In the light of the previous results, we tried to account for the 
observational biases involved in measuring z that come into 
play when one considers the samples of unknown-z GRBs with 
some evidence for anisotropy, S2 and SnEnZ. To this aim, us¬ 
ing method 1 we generated another set of 1000 synthetic sam¬ 
ples with the same size as S2, which took into account two main 
sources of bias: 
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Fig. 4. Fraction of GRBs with unknown redshift as a function of Galac¬ 
tic dust extinction. 

- Galactic dust extinction disfavours optical observations of 
low-Galactic-latitude GRBs; 

- ground-based follow-up observations are less probable 
around the celestial poles (iFvnbo et al.ll2009t) . 

To model the first bias, we considered the G alactic extinction as 
measured bv ISchlaflv & Finkbeined (1201 Ih in terms of AvQ We 
then obtained the Av(Gal) distribution for the positions of all the 
GRBs belonging to STot and split it into two groups, with and 
without measured z respectively. For each Av(Gal) bin we cal¬ 
culated the fraction of GRBs with unknown z over the total. The 
result is shown in Fig. |4l As long as it is Av(Gal) < 1 mag the 
fraction is about 2/3, i.e. the same value that is obtained over the 
entire Swift catalogue. However, for Av(Gal) > 1 mag almost all 
of the GRBs have no redshift. In generating the synthetic sam¬ 
ples analogous to S2, each random position was accepted de¬ 
pending on the outcome of a binomial variate: the acceptance 
probability was set to 1 if Ai/(Gal) > 1 mag, to 2/3 otherwise. 

To model the effect of the second bias, we preliminarily 
compiled a list of the ground facilities that mostly contributed 
to measuring GRB redshifts. For each declination, we consid¬ 
ered the visibility from each location by calculating the fraction 
of time that a generic source at the given declination spends at 
> -1-30° above the local horizon. Such threshold value matches 
average observational constraints. We then averaged it over the 
full set of locations, weighting by the multiplicity of the num¬ 
ber of telescopes at a given location (e.g., in La Palma there are 
several telescopes which provided GRB redshift measurements). 
We have therefore come up with a probability of having z mea¬ 
sured as a function of declination. This bias was incorporated in 
the MC simulations as a further criterion for a random position 
to be accepted, depending on the outcome of a binomial variate: 
the acceptance probability was set to the complement to one of 
the probability of having z measured for that declination. 

We also applied method 2 by weighting each sky pixel by the 
product of the sky exposure and the two above probabilities. We 
then calculated the mean value for each moment and the corre¬ 
sponding uncertainties as in Section lT^ under the same assump¬ 
tions. The results are reported in Table |2l Both methods gave 
fully consistent results as was the case for the other samples. 

It is noteworthy that the deviation from isotropy of the Galac¬ 
tic quadrupole moment (sin^ b - 1/3) tapers off from -2.8cr to 
-1.8cr when observational biases are accounted for. Likewise, 


the equatorial quadrupole moment (sin^ <5-1 /3) decreases from 
- 1 - 2 .3cr to -i-L7cr. 

We also considered the Galactic dust correction alone and 
found that it accounts for -90% of the total change in the ex¬ 
pected Galactic quadrupole moment, so it dominates over the 
correction for the locations of ground telescopes. 

4.5. Constraining the fraction of apossibie Gaiactic 
popuiation 

The Galactic quadrupole moment (sin^h - 1/3) of unknown- 
z events is marginally deviating from isotropy: this holds both 
for the FREDs (S2) and, to somewhat lower extent for the non- 
FREDs (SnEnZ). This would imply an excess of events at low 
Galactic latitudes. Taken at face value, when one overlooks the 
impact of the observational biases discussed in Section l4Al and 
plainly assumes that this is mostly due to the presence of a 
Galactic population as T12 did, one can estimate its fraction 
within the sample of unknown-z FREDs. 

We did so by generating 100 synthetic samples with as 
many events as in S2, a fraction of which was obtained from 
an isotropic distribution, and the remaining one from a random 
selection of a sam ple of 334 known Galac tic sources from the 
BAT catalogue by ICusumano et al.l (l2010l) PI Both populations 
were selected compatibly with the exposure map. 

By varying every time the fraction of the Galactic compo¬ 
nent, we ended up with the best match between observations and 
simulations. This turned into an estimated fraction of Galactic 
intruders in S2 of / = (27 + 10)% (Icr), similar to T12’s conclu¬ 
sions. In practice, given that the evidence for anisotropy is only 
marginal, we provided a 90%-conhdence upper limit of / < 40% 
by demanding that at least 90% of the synthetic samples had a 
non-zero quadrupole moment at >3cr. 

However, neglecting the observational biases has the effect 
of overestimating the weight of the possible Galactic population. 
We therefore repeated the same test by accounting for those bi¬ 
ases. Eollowing the procedure of Section [AT] we came up with 
a more realistic estimate of / = (19 + 11)% with a 90% upper 
limit of / < 34%. 

5. Discussion 

We found no evidence for a signihcant dipole moment in the 
different FRED samples observed by Swift. However, we point 
out that this does not clash with T12’s claim: in contrast with 
what they assert, cos b is by no means a dipole, since a dipole 
moment is the cosine of the angle between a generic direction 
and a fixed one in the sky, not a plane. Thus, the tension is just 
in the claim but not in the results, which do agree. In particular, 
signihcant dipole moment was found neither along the Galactic 
poles (sin b), nor towards the Galactic centre (cos b cos /)• 

In contrast, and in agreement with T12, we found that the 
Galactic quadrupole moment (sin^ h - 1 /3) of unknown-z events 
marginally deviates from isotropy, implying an excess of events 
at low Galactic latitudes. 

However, before giving credence to the possible existence of 
such a relatively abundant disguised Galactic community, one 
must carefully evaluate the impact of observational biases other 
than the non-uniform exposure map. That not just FREDs, but 


^ We made use of known sources just to mimic in a credible way 
what the distribution of a Galactic population should look like. Clearly, 
such possible unrecognised Galactic sources must obviously be differ- 
Available at http://irsa. ipac. caltech. edu/applications/DUS’Efi : from the known ones. 
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all GRBs with unknown z tend to cluster towards the Galactic 
plane, likely suggests that the impact of the process of mea¬ 
suring z might be far from negligible. This possibility appears 
to be corroborated by the absence of non-zero quadrupole mo¬ 
ments in the entire sample of BAT events (STot), in spite of 
its having the best statistical sensitivity. Low-Galactic-latitude 
GRBs are notoriously disfavoured in their optical counterpart 
detection and consequent z measurement because of the strong 
Galactic dust extinction. Moreover, our result that unknown-z 
events tend to cluster around the celestial poles (|d| > 60°), 
where ground-bas ed follow-up observations are less probable 
(iFvnbo et alj|2009l) . lends further support to the z-related obser¬ 
vational bias explanation. In the light of these considerations, 
we modelled the impact of both effects in shaping the unknown- 
z sample S2 by preliminarily studying how the measurement of z 
is hampered by the amount of Galactic dust content for the entire 
sample of Swift GRBs. Furthermore, we modelled the probabil¬ 
ity for a generic source at a given declination to be followed 
up spectroscopically by estimating its visibility from the most 
frequently used ground facilities. As a result, we found that the 
deviation from isotropy of the Galactic quadrupole moment de¬ 
creased from -2.8cr to -1.8cr, thus deflating the body of evi¬ 
dence of a significant contamination due to a disguised popula¬ 
tion of Galactic transients. 

We delved deeper into such matter by estimating the con¬ 
taminating fraction of a Galactic population that would yield the 
same Galactic quadrupole moment as the one observed in the 
sample of unknown-z FREDs. When the observational biases are 
properly accounted for, we estimate / = (19+ll)%(lcr) with a 
90% upper limit of / < 34%. 

To gain further insight, we examined the spectral and tempo¬ 
ral properties of the y-ray prompt emission along with the X-ray 
afterglow of the S2 sample as a function of the Galactic position 
and compared their distributions with those derived for a sample 
of cosmological GRBs with measured z. Specifically, we esti¬ 
mated the rise and the decay times of the y-ray profiles as well 
as the photon index in the 15-150 keV range of the spectrum 
extracted over the Tgo interval. The X-ray afterglow properties 
were tak en from the classif i cations provided by S wift-XRT cat¬ 
alogues (lEvans et al.l l200^lMargutti et al.ll2013h and from the 
online Leicester cataloguqj for the recent GRBs. Eor none of 
these observables we found evidence for a different distribution 
between S2 and cosmological GRBs. 

Noteworthy is that the fraction of short-duration GRBs in 
S2 is less than 7% (8/119), i.e. far from enough to account for 
the deviation from isotropy that we observe. Several investi¬ 
gations on past short an d very short GRB catalogues brought 


evidence for anisotronv dBalazs et al 

Il998t Magliocchetti et al.l 

120031:1 Vavrek et al.ll2008l; ICline et al.1 

2003112005b. However, the 


low fraction of short GRBs in our Swift samples rules them out 
as the main cause for the marginally observed anisotropy. 


6. Conclusions 

Triggered by previous results, we searched for the possible pres¬ 
ence of disguised Galactic transients in the catalogue of GRBs 
detected by Swift-BAfT up to September 2014. We focused on 
ERED GRBs, given that their light curves are often observed in 
a broad range of Galactic high-energy transients. The search was 
based on the statistical analysis of the sky distribution both with 
reference to our Galaxy and in a reference-system-independent 
way. While no deviation from isotropy was found in the dipole 

^ http://www.swift.ac.uk/ 


moment, we found marginal (2.8cr) evidence for a non-zero 
Galactic quadrupole moment for EREDs with unknown redshift 
z, suggesting a clustering at low Galactic latitudes, in partic¬ 
ular towards the celestial poles (|(i| > 60°). Should this pos¬ 
sible clustering on the Galactic plane be due the presence of 
an unidentified Galactic population, we constrain its fraction to 
/ = (27 + 10)%. 

However, non-ERED GRBs with unknown z also show a 
similar behaviour, even though with lower statistical significance 
(2. Icr). If one drops the unknown-z condition, such evidence dis¬ 
appears. It is known that there are observational biases that dis¬ 
favour the measurement of z: the large dust extinction towards 
the Galactic plane along with the lower probability of ground- 
based observations of polar (|d| > 70°) events, can contribute 
to enhance the degree of anisotropy. Once such biases are ac¬ 
counted for, the evidence for a clustering near the Galactic plane 
decreases from ~ 3cr to ~ 2cr, thus weakening the claim for a 
Galactic contamination. Eurthermore, comparing the properties 
of the y-ray prompt emission as well as of the X-ray afterglow 
with those of GRBs with measured z did not reveal any statisti¬ 
cal difference. Hence, our results suggest marginal evidence for 
the presence of an unrecognised Galactic population among the 
GRBs with unknown redshift. 

In conclusion, when z-measurement-related biases are ac¬ 
counted for, the fraction of a possible disguised Galactic pop¬ 
ulation of high-energy transients among Swift-BAfT GRBs with 
unknown redshift is estimated / = (19 + 11)%, with a 90%- 
confidence upper limit of / < 34%. 
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